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Abstract 

Melanocytic nevi are benign proliferations that sometimes turn into malignant melanoma in a way that is still unclear from 
the biochemical and genetic point of view. Diagnostic and prognostic tools are then mostly based on dermoscopic 
examination and morphological analysis of histological tissues. To investigate the role of mechanics and geometry in the 
morpholgical dynamics of melanocytic nevi, we study a computation model for cell proliferation in a layered non-linear 
elastic tissue. Numerical simulations suggest that the morphology of the nevus is correlated to the initial location of the 
proliferating cell starting the growth process and to the mechanical properties of the tissue. Our results also support that 
melanocytes are subject to compressive stresses that fluctuate widely in the nevus and depend on the growth stage. 
Numerical simulations of cells in the epidermis releasing matrix metalloproteinases display an accelerated invasion of the 
dermis by destroying the basal membrane. Moreover, we suggest experimentally that osmotic stress and collagen inhibit 
growth in primary melanoma cells while the effect is much weaker in metastatic cells. Knowing that morphological features 
of nevi might also reflect geometry and mechanics rather than malignancy could be relevant for diagnostic purposes. 
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Introduction 

Melanocytic nevi are benign proliferations of melanocytes, the 
skin cells that produce the pigment melanin. They are by 
definition benign, but 50% of malignant melanomas arise from 
pre-existing nevi. Current diagnostic methods of melanocytic 
lesions are based on hystopathology and dermoscopy which reveal 
wide morphological diversity and evolution patterns whose origin 
is still unknown [1,2]. Leaving aside the controversial cases of Spitz 
and blue nevi, nevi are categorized from dermoscopic analysis into 
globular, reticular, structureless brown and mixed patterns (the 
latter are subdivided into mixed pattern with globular structureless 
areas in the center or mixed pattern with globuls at the periphery) 
[3]. The classical theory of nevi evolution formulated in 1893 by 
Unna [4] claims that nevi originate from melanocytes proliferating 
at the dermo-epidermal junction (lentigo simplex and junctional 
nevus) that form nests (compound nevus) and eventually migrate 
completely into the dermis (dermal nevus). More recently, Cramer 
described an opposite model known as the theory of upward 
migration [5] . Cramer suggested that precursor cells of melano- 
cytes deriving from pluripotent stem cells of the neural crest 
wander during embryogenesis along nerves into the dermis, 
mature here and finally migrate as functional melanocytes into the 
epidermis. Cramer's last migration step, however, was never 
broadly accepted. More recently Kittler et al. [6] suggest that nevi 
can migrate not only vertically but also horizontally, explaining 



how nevi can expand in the course of time. The underlying 
limitation of these theories is that they are based on hystopatho- 
logical observation only, and do not reflect the dynamics of an 
individual nevus. From a diagnostic point of view it would be 
extremely useful to correlate the morphological features of nevi to 
the degree of malignancy or pre-malignancy. This issue is intricate 
because nevi proliferate in a complex microenvironment that can 
mediate cell behavior through the composition, structure, and 
dimensionality of the extracellular matrix (ECM), the polymeric 
scaffold that surrounds cells within tissues. For example, a 
malignant phenotype can be reverted into a nonmalignant one 
by specifically blocking aberrant adhesion of the cancer cell to its 
extracellular scaffold [7]. 

Recent research shows that mechanical properties of the tumor 
microenvironment and of the normal tissue can influence tumor 
growth and dynamics, in way that is still poorly understood [8— 
1 3] . Mechanical stresses such as those experienced by cancer cells 
during the expansion of the tumor against the stromal tissue have 
been shown to release and activate growth factors involved in the 
progression of cancer [14]. Moreover, the stiffness of the matrix 
surrounding a tumor determines how cancer cells polarize, adhere, 
contract, and migrate, and thus regulates their invasiveness [15]. 
Forces exerted by cancer cells as they migrate through the ECM 
has been quantified accurately using traction force microscopy 
[16,17] Another possibility is that mechanical stresses directly 
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regulate the growth and death rates of cancer cells as shown by 
Montel et al. who induce osmotic stress by adding dextran, a 
biocompatible polymer that is not metabolized by cells [11,12]. 
Several studies in the literature indicate important changes in 
cellular functioning due to osmotic pressure [18,19], but the 
stresses involved (in the MPa range) were orders of magnitude 
larger than those (in the kPa range) studied in Refs. [11,12]. It is 
interesting to notice that compressive stresses of slightly less than 
1 kPa applied through a piston were recently found to induce a 
metastatic phenotype in cancer cells [13]. While osmotic pressure 
may have a different origin than compressive mechanical pressure 
the effect on cells is exactly the same. This was demonstrated 
experimentally in Refs. [11,12] by applying osmotic pressure 
directly on the cells, by adding dextran to the growth medium, or 
indirectly, by placing the cells inside a dialysis bag which was then 
placed into a dextran solution. In both cases, the effect of pressure 
on the cells was found to be exactly the same. This is not 
surprising, since osmotic pressure corresponds indeed to a real 
mechanical pressure on the membrane due to the collisions with 
solute molecules. 

In a recent paper, Simonsen and coauthors showed that 
interstitial fluid pressure (IFP) was associated with high geometric 
resistance to blood flow caused by tumor-line specific vascular 
abnormalities in xenografted tumors from two human melanoma 
lines with different angiogenic profiles [20]. In another recent 
paper, Wu and colleagues investigated how nonlinear interactions 
among the vascular and lymphatic networks and proliferating 
tumor cells might influence IFP transport of oxygen, and tumor 
progression [2 1] . They also investigated the possible consFuences 
of tumor-associated pathologies such as elevated vascular hydrau- 
lic conductivities and decreased osmotic pressure differences. All 
these parameters might affect microenvironmental transport 
barriers, and the tumor invasive and metastatic potential, opening 
interesting new therapeutic approaches. In general, understanding 
the influence of mechanical stress on cancer growth could shed 
new light on tumor development and progression. 

To investigate the role of the environment in the growth of 
melanocytic nevi, we introduce a computational model of elastic 
cells proliferating in a layered non-linear elastic medium 
representing the skin. Several computational models have been 
introduced in the past to simulate cancer growth in vitro and in 
vivo. They range from individual cell models [22-26], to lattice 
models [27-29] and continuum models [29-31]. The aim of our 
investigation is to correlate the morphological features of the 
growing nevus with the location of the originating melanocyte and 
the mechanical properties of the environment. To this end, a 
natural strategy is provided by individual cell models where each 
cell can be mechanically deformed and moves according to 
mechanical forces due to other cells and to the ECM. Models of 
this kind have been used in the past to simulate the growth of the 
epidermis [32,33] and the growth of cancer cell colonies inside a 
confining medium [26]. 

We first report experimental results on melanoma cell lines 
showing that a low osmotic pressure reduces proliferation, in 
agreement with other reports for other tumors [11,12]. We have 
performed experiments on two cell lines, the first obtained from a 
primary tumor and the second from a metastatic one, both from 
the same patient. We find that melanoma cells stemming from a 
primary tumor are less sensitive to pressure than the metastatic 
ones. We find a similar results by growing the same cells in 
collagen coated plates, simulating the effects of the ECM. 
Melanoma is widely believed to be originating from melanocytes, 
which are not easily cultivated in vitro. Therefore melanoma cell 
lines is the most sensible in vitro model for nevi. We expect that 



the effect of pressure on proliferation seen in melanoma cells will 
be even stronger for melanocytes since melanoma cells are much 
more resistant. 

Numerical simulations of the model yield different growth 
patterns for nevi originating in the epidermis, either close to the 
basal membrane or to the stratum corneum, or in the dermis. 
Different initial locations give also rise to different mechanical 
properties as we show by measuring compressive stresses on the 
proliferating melanocytes. Compressive stresses are fluctuating 
strongly inside the nevus, with an average stress that typically 
increases until the ECM breaks leading to stress drop. A crucial 
role in the growth pattern is played by the basal membrane, 
separating the epidermis from the dermis. Depending on its 
mechanical properties, the membrane can resist the pressure 
induced by the nevus or fail leading to the invasion of the dermis. 
In this case, invasion is not correlated with an increased 
malignancy of the tumor but just to mechanics. This observation 
of course does not exclude that malignant phenotypes may be 
induced by mechanical stresses as previously reported in the 
literature [13]. It is in fact known that malignant cells break the 
basal membrane by matrix metalloproteinases (MMPs). We 
simulate this process explicitly showing that it accelerates the 
invasion of the dermis, but the general morphology remains the 
same. 

Model 

We simulate an individual cell based model [24-26] for the 
growth of melanocytic nevi in the skin. We consider a two 
dimensional vertical section of the skin (see Figure 1 a) consisting of 
two main layers, the epidermis and the dermis, separated by the basal 
membrane. The deepest layer of the skin, the hypodermis, is not 
considered in the present study. In our model, the dermis 
constitutes the bottom part of the skin, whereas the upper layer 
is represented by the epidermis and by its most superficial part, the 
stratum corneum. The model describes the skin as a highly organized 
elastic organ, whose mechanical properties are specific of any of its 
constituent layers. Periodic boundary conditions are applied to any 
of the skin layers along the horizontal direction, while open 
boundary conditions are imposed in the vertical directions. 

Dermis 

The dermis is arbitrarily divided into two anatomical regions: 
the papillary and reticular dermis. We consider only the papillary 
layer, whose thickness we set to 150/mi. We model the dermis as a 
network of cells with an average radius ~15/mi, arranged in a 
disordered lattice and connected by non-linear springs. The 
disordered lattice is obtained starting from a triangular lattice, 
displacing the nodes randomly by a small amount. The resulting 
configuration is triangulated by a Voronoi algorithm. Nearest 
neighbors cells are connected by non-linear springs with a rest 
length equal to the resulting distance between the displaced cells. 
Cells here represents fibrolasts, macrophages, and adynocytes 
present inside the papillary layer, while springs model the effect of 
the ECM, composed by collagen fibers and elastin. We set the 
average length of these springs at ~60/mi. The dermis also 
contains several irregularities such as blood vessels, lymph vessels, 
nerve endings and skin appendages such as air follicules, small hair 
muscles and sebaceous glands, but we ignore these details in the 
present model. The mechanical behaviour of the dermis has been 
studied experimentally using in-vivo suction and tension tests as 
reported in [34-36]. In these studies, the mechanical response of 
the dermis layer to an external stress was found to obey to the 
following non-linear James-Green-Simpson pressure-displacement 
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Figure 1. Computational model, a) A two dimensional section of the skin (left image courtesy of Dr. Claudio Clemente) is simulated as a non- 
linear elastic layered material (right). The top layer is the stratum corneum (pink cells), resting on the epidermis (blue cell). The epidermis is separated 
from the dermis (green cells) by a basal membrane (grey). In a typical simulation we consider a melanocyte (yellow) dupicating inside the skin, either 
in the epidermis, or in the dermis, b) The connecting fiber have a non-linear elastic behavior with a fracture strength set at a Y = 1.8kPa. c) When in 
contact, the cells interact by a finite-thickness Hertz law. The corresponding zero-thickness law is reported for comparison. 
doi:10.1371/journal.pone.0094229.g001 



relationship [37] 

a = 2C w (a- 1) +6C„ (k 2 -k- 1 + 1 + 1 - 1) (1) 

where o and X represent respectively the stress and the stretch of 
the tissue, the strain being just £ = 2—1. The expression (1) is 
derived from the strain energy density function usually adopted to 
mimic the non-linear behaviour of the skin [38] . Here we use the 
same law to model the deformation of the ECM springs, setting 
their elastic constant to C\q =0.1 125kPa and C\\ =0.315kPa, 
and introduce a failure stress at o y = 1.8kPa, at which the spring 
breaks. The non-linear elastic behavior of the springs is illustrated 
in Fig lb, showing that when the stress overcomes Oy the spring 
breaks and the stress goes to zero. 

Epidermis 

The epidermis is a multilayered tissue composed of the stratum 
granulosum, stratum spinosum and stratum corneum. Since the 
elastic properties of the stratum corneum are cosiderably different 



from the other two, we will consider it separately. The thickness of 
the epidermis (stratum granulosum and stratum spinosum) is 
150/xm and it is composed mainly by keratinocytes, whose average 
radius is ~20/mi. In the model, keratinocytes are placed at the 
nodes of a disordered lattice, obtained as discussed above for the 
dermis. The turnover rate of melanocytes is assumed to be much 
faster than the natural turnover rate of keratinocytes, modelled 
explicitly in Ref. [32] but ignored here. This is justified since the 
typical turnover for melanoma cells is 1-2 days [42] which is much 
smaller than the turnover of keratinocytes (50-70 days). To 
simulate the packing of keratinocytes in the epidermis, we chose 
the radius of each cell after the Voronoi triangulation so that no 
neighboring cell overlaps. Cells are then coupled by non-linear 
springs accounting for the neo-Hookean behavior observed in in- 
vivo experiments [34-36] . Experiments suggest that the non-linear 
pressure-displacement relation in Eq. (1) captures the mechanical 
properties of the epidermal layer with a good accuracy using the 
same parameters values obtained for the dermis [34,35]. In the 
model, we link the cells by non-linear springs with 
Ci 0 = 0.1125kPa, Cn=0.315kPa and a r = 1.8kPa. The average 
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length of the links between cells coincides with the average cell 
radius, hence two keratinocites may often be in contact and will 
thus interact as we discuss below. 

Stratum Corneum 

The stratum corneum is composed of dead corneocytes of an 
average radius of ~20/mi. Although it may consists of up to 15—20 
layers of corneocytes, in our simplified numerical model only one 
single layer is taken into consideration. Experiments [34-36] have 
provided evidence that the stratum corneum mechanical proper- 
ties are consistent with a linear Hookean law. Here, we consider 
elastic cells linked by linear springs with Young's modulus 
E=\.5kPa and failure stress set to cry = 2kPa. The stratum 
corneum cells are linked to the epidermis by linear springs whose 
Young's modulus is E= 1.6125kPa. 

Basal membrane 

The basal membrane constitutes and elastic sheet separating the 
dermis from the epidermis. It has a complex structure composed of 
the lamina reticula and lamina basale connected by collagen 
fibrils; usually the thickness is around 0.5— In our model the 
basal membrane is represented by an assembly of fictitious cells of 
radius 0.5/mi, connected by linear springs with Young's modulus 
which we set either 0.15kPa or 0.3kPa. The failure stress is of the 
springs is set to Oy = 2kPa. The fictitious cells composing the basal 
membrane have Hertzian repulsive interactions with the other 
cells thus providing an effective barrier to their motion. To 
represent the highly corrugated structure of the membrane, we use 
a simple periodic function with amplitude equal to 1.5/im and 
period 30/mi. Dermis and epidermis cells are connected to the 
basal membrane from below and above respectively, these springs 
obey to a linear Hookean force with Young's modulus equal to 
£' = 0.18725kPa. In the model, the basal membrane can be 
broken mechanically, but we also consider the production of 
MMPs. To model this effect, we stipulate that at each time step 
each cell composing the nevus can dissolve one link with a 
probability that decreases with the distance from the cell: 



P MMP =Pbreak 



with pbreak = 10 2 and r break = 30/mi. 



^break 



(2) 



Cell mechanics 

Cells are modelled as non-linear elastic spheres and, for the sake 
of simplicity, we employ the same mechanical constants for all the 
cells composing the epidermis the dermis and the nevus. When 
two cells are in contact they repel by a finite-thickness Hertzian 
law [39] 

F Hertz = 3 (1—2) VIM , (3) 

where E is the cell's Young's modulus, set to E~ lkPa [40], v is 
the Poisson's ratio, R = 



R\R2 
Rl+Rl 



with R\ and R2 radii of the 



contact cells, 5 is the indentation depth and co = 0.1 — (see Fig. lc). 

We notice that cells are often modelled as incompressible (v = 0.5) 
[40], but the value of Poisson ratio is difficult to estimate 
experimentally. Individual cell models use compressible cells as we 
do and employ v<0.5 [26,41], Here we use v = 0.33 but its precise 
value is not really important in the framework of our simulations: 



changing v would just imply a modification of the pre-factor of Eq. 
3, leaving qualitatively unaffected the numerical results hereby 
displayed. We notice that, although Eq.(3) refers to the experi- 
mental situation of a microsphere indenting a finite-sized soft gel, 
it captures the non-linear corrections to the usual Hertzian law, 
needed when large deformations occur. In most simulations, nevi 
cells are interacting by Hertzian repulsion only, while the other 
cells are typically linked also by non-linear spring. We have also 
tested for the effect of adhesion between nevi cells introducing a 
radial force /adhesion = kadd, with k a d = 1.5 x 10 -5 N/m. Adhesive 
forces only act when the indentation depth 5 is larger than 
— 20/mi. 

Interactions with the ECM 

In the simpler implementation of the model, nevi cells have no 
direct interaction with the ECM. An indirect interaction exist since 
nevi cell interact with other cells which are held together by the 
ECM. This indirect interaction is, however, weak inside the 
dermis. To overcome this limitation, we also consider a direct 
interaction between nevi cells and the ECM. In practice, we follow 
the same strategy used to model the basal membrane and place 
three fictitious cells with radius 5^111 on each ECM link. In this 
way, nevi cells interact directly with the ECM. 

Growth dynamics 

The quasi-static growth of the nevus is simulated by randomly 
selecting a single cell among those belonging to the dermis or 
epidermis. This cell represents the first proliferating cell and 
therefore does not present any linkage with the surrounding skin 
network. The growth process is obtained by duplication of the 
initial cell, so that the new cell is still in contact with the first one, 
but its growth direction is picked at random, i.e. lattice-free: such a 
choice is motivated by the fact a proliferative population will 
always eventually fill the lattice, suppressing cell-to-cell crowding 
effects [29] . The insertion of a new cell encompasses an overall 
rearrangement of the entire network, which is achieved by 
minimization of the stress within the system. Since growing cells 
are completely disconnected from the surrounding skin tissue, they 
only interact through Hertzian-like forces (Eq. 3) and possibly 
through adhesion. Two different choices of protocol are applied to 
the duplication mechanism: a cell can either be picked randomly 
among those composing the nevus, or the cell can be selected with 
a rate that depends on the compressive stress. In particular, the 

probability Pf(N) cce « N » is assigned to each cell, where (?i(N) 
represents compressive stress of the z-th cell and <cr(7V)> is the 
average stress. Finally, the cell that duplicates is selected with 
probability Pi(N). 

Results 

Experiments 

Osmotic pressure affects the growth of primary 
melanoma cells. We test the effect of low osmotic pressure 
(from 0. 1 kPa to 1 kPa) on cellular proliferation of primary and 
metastatic human melanoma cell lines, IgR39 and IgR37, 
respectively at short (3 days) and long time (6d ays). As shown in 
Fig. 2, the treatment with 1 kPa for 6 days decreases the cellular 
proliferation of IgR39 while IgR37 are only slightly affected. 
According to Table 1 , the osmotic pressure increases cell death to 
necrosis while the percentage of cells in early apoptosis slightly 
decreases at low pressures and then returns to the original value at 
higher pressure (Table 1). We are tempted to attribute this last 
result to a statistical fluctuation. 
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Figure 2. Effect of osmotic pressure on cell proliferation in melanoma. 500 cells/well were plated on 96 multiwells (lgR39 and lgR37). The 
day after plating the cells were submitted to different osmotic pressure (from 0.2 to 1 kPa) for 3 or 6 days. At the end of the incubation the cells were 
fixed with 50%TCA for 2 hours at 4C and air-dry at room temperature. Thus, the cells were incubated with 0.05% SBR solution for 30 minutes at room 
temperature and then quickly rinsed four times with 1% acidic acid to remove unbound dye. Finally, the protein-bound dye was solubized with 
10 mM TRIS and OD was measured at 510 nm with microplate reader 550 (Bio-Rad). a) The growth of lgR39 (non-metastatic) cells is not affected by 
pressure after 3 days, b) but cells grow significantly less after 6 days, c) lgR37 (metastatic) cells are unaffected by pressure both after c) 3 or d) 6 days. 
Statistically significant results according to the KS test (/?<10~ 2 ) are denoted with * 
doi:1 0.1 371 /journal. pone.0094229.g002 



Table 1. Effect of osmotic pressure on cell death. 
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necrosis or early apoptosis after 3 or 6 days for various osmotic pressures in a typical experiment. 
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We also perform a colony formation experiment using the 
crystal violet assay on primary and metastatic melanoma cells, 
following the prescriptions of Ref. [42] . In Fig 3a we report the 
cumulative colony size distribution P(s), where s is the number of 
cells in each colony, for IgR39 cells with or without osmotic 
pressure. Each distribution is obtained by combining the result of 
colonies obtained in six different wells. Typical images of the 
colonies are reported in Fig 3c. We fit the distributions with the 
solution of a continuum time branching process in which each cell 
divides with rate y which yields a cumulative colony size 
distribution after t days given by [42] 

P(s) = {\-e- yt ) s for s>\. (4) 

From the fit, we estimate the average division rate as y = 0.60 
divisions /day in normal conditions and y = 0.51 divisions /day 
under an osmotic pressure of 1 kPa. In Fig. 3b we report the 
average colony size with its standard error for IgR39 and IgR37 
cells with and without osmotic pressure. While osmotic pressure 
hinders colony growth in both cases, the effect is much stronger for 
primary tumor cells than in metastatic ones. 

Effect of collagen on the growth of melanoma cells. To 
confirm that the general effect of stress we observe on growing 
melanoma cells is not restricted to osmotic pressure, we perform 
proliferation experiments using collagen coated plates simulating 
the effect of to the extracellular matrix on the cells. As shown in 
Fig. 4, the presence of collagen significantly reduces proliferation 



in primary melanoma cells (IgR39) while the effect is much weaker 
in metastatic cells (IgR37). We next induce an osmotic pressure of 
1 kPa in the collagen coated plates and find no additional 
significant result with respect to the case with collagen alone. We 
can interpret this result by assuming that collagen already induces 
a strong stress on the cells so that additional osmotic stress has no 
significant effect of the growth. We can thus suggest that osmotic 
pressure and collagen both act in similar ways on melanoma cells: 
they reduce proliferation but the effect is much stronger in primary 
cells than in metastatic ones. 

Numerical simulations 

Morphology of melanocytic nevi depends on the location 
of the initiating cell. We perform numerical simulations of the 
growth of melanocytic nevi for different conditions. In particular, 
we study the morphology resulting from the growth using different 
locations of the initial seed. As illustrated in Fig 5, we consider four 
type of initial locations: deep inside the epidermis (Fig. 5a and 
movie SI), in the epidermis but close to the basal membrane, 
either in the minima (Fig. 5b, Fig. 5e and and movies S2-S3) or on 
the maxima of the membrane (Fig. 5c, Fig. 5f and and movie S4- 
S5), or inside the dermis (Fig. 5d and movie S6). For initiating cells 
close to the basal membrane, we consider two different scenarios 
with a strong (Fig. 5b and 5c) or weak membrane (Fig. 5e and 5f). 
Fig. 5 indicates that the morphology of the nevus depends 
substantially on the location of the initiating cell. When the nevus 
starts up in the epidermis the growth is mostly horizontal and far 
from the dermis (Fig. 5a). For initiating cells residing close to a 




s [number of cells] 

Figure 3. Effect of osmotic pressure on colony formation in melanoma, a) The cumulative distributions of colony size obtained from lgR39 
cells under 1 kPa osmotic pressure with respect to the control (0 kPa). The curves are the fit with a continuous time branching process model (see Ref 
[42]) yielding a division of rate of y = 0.60 that is reduced to y = 0.51 under osmotic pressure, b) The average value of the colony size distribution with 
the associated standard error for lgR39 and lgR37 cells. Statistically significant results according to the KS test (/?< 10~ 2 ) are denoted with *. c) The 
images show two representative examples of the colonies for 0 kPa and 1 kPa conditions. 
doi:10.1371/journal.pone.0094229.g003 



PLOS ONE | www.plosone.org 



6 



April 2014 | Volume 9 | Issue 4 | e94229 



Mechanical Properties of Growing Melanocyte Nevi 



lgR39 



lgR37 



2 - 



E 
c 
o 



Q 
O 



1 - 



* 

53% 



* 

18% 



~ i i i r~ 



i i i i 



~~ i i i r~ 



~i i i r~ 



~i i i r~ 



~i i i r~ 



OkPa +Collagen +Collagen OkPa 
+ 1 kPa 



♦Collagen +Collagen 
+ lkPa 



Figure 4. Effect of collagen coating on cell proliferation in melanoma. 500 cells/well were plated on 96 collagen coated multiwells (lgR39 
and lgR37) and submitted to 1 kPa osmotic pressure for 6 days. Cell growth was evaluated as in Fig.2. The graph shows that the presence of collagen 
slows cell growth both in primary (lgR39) and metastatic (lgR37) cells but the effect is much stronger for primary cells. Additional osmotic pressure 
does not significantly alter the results. Statistically significant results, according to the KS test (p<10~ 2 ), are denoted with * 
doi:10.1371/journal.pone.0094229.g004 



strong basal membrane the growth follows the profile of the 
membrane (Fig. 5b) or expands vertically in the epidermis (Fig. 5c) 
depending on whether the initiating cells lies in the minima or 
maxima defined by the membrane profile. When the membrane is 
weak, however, it is easily broken by the growing cells that then 



invade the dermis (Fig. 5e and Fig. 5f). The growth is more radial 
when the initiating cells is in the minimum of the basal membrane 
(Fig. 5e) and more vertical when it lies close to the maximum 
(Fig. 5e). 




Figure 5. Morphology of nevi for different locations of the initiating cell for random growth. Illustration of the results of numerical 
simulations for nevi grown from melanocytes located in different positions in the skin and for different mechanical properties of the basal membrane. 
Here growth occurs randomly and independently on compressive stress. Growing melanocytes are shown with a varying color that reflects their 
compressive stresses o according to the color bar. The maximum a max and minimum c- mi n for each configuration are equal to (in kPa): a) 0.61 -0.77, 
b) 0.0043-0.21, c) 0.085-0.57, d) 0.0013-0.149, e) 0.0047-0.37, f) 0.0067-0.72. a) Nevi grown from melanocytes residing inside the epidermis 
tend to grow horizontally and do not spread towards the basal membrane, b) Nevi grown on the minima of a strong basal membrane tend to grow 
roughly parallel to the membrane itself, c) Nevi growing from the maxima of a strong basal membrane tend to grow vertically in the epidermis, d) 
Dermal nevi tend to have radial shape. When the basal membrane is weak, e) nevi growing from the minima of the basal membrane invade the 
dermis in a radial fashion, while f) when they start from maxima of the basal membrane the invasion of the dermis occurs more vertically. 
doi:10.1371/journal.pone.0094229.g005 
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Mechanical stresses in growing melanocytic nevi. We 

record the compressive stresses sustained by the melanocytes 
during the growth of the simulated nevi. A pictorial view of the 
stress in each cell is shown in Fig. 5, where the color scale depicts 
compressed cells with a scale going from yellow (low stress) to red 
(high stress). The images show that stresses are in general 
heterogeneously distributed and, as illustrated by the supplemen- 
tary movies, fluctuate considerably after each division step. We 
notice that typically the most stressed cells are located in the inner 
part of the nevus, reproducing a general feature observed in tumor 
spheroids [11,12]. Furthermore, we see the formation of stress 
chains, typical of granular assemblies [43] . To better quantify the 
evolution of the compressive stress we measure the average 
compressive stress a experienced by the cells as a function of the 
number of cells in the nevus. The results reported in Fig. 6 
correspond to the different cases explored in Fig. 5. In all cases, 
pressure builds up as cells are duplicating until at some point the 
connective tissue breaks drastically reducing the average compres- 
sive stress experienced by the cells. This effect is particularly clear 
in the case of a nevus growing in the epidermis (see Fig. 6a). The 
large error bars reported in Fig. 6 indicate that stresses vary widely 
between individual cells but also among different realizations of 
the process. 

To better quantify the stress fluctuations experienced by 
individual cells, we report in Fig. 7 an estimate of the distribution 
(probability density function) p(cr) of the compressive stresses <7 for 
different initial conditions. The distribution is sampled over 
different realizations and over different cells. Since the distribution 
is expected to change with the number of cells present in the nevus 
(as shown from the error bars in Fig. 6), we consider only the case 
for which the average pressure is highest (the peak stress in each 



curve in Fig. 6). In all cases the distribution displays an 
approximately exponential tail (i.e. p(<j)— exp ( — o/gq)). 

Compressive stresses are much higher for nevi growing in the 
epidermis than in the dermis. This is due to the general 
assumption made in the model that cells are more packed in the 
epidermis than in the dermis. When we include a direct interaction 
between the nevi cells and the ECM, as discussed in the model 
section, compressive stresses in the dermis increases but the 
general features of the results are unchanged (see Fig. 8 and Movie 
S7). We have also tested for the effect of adhesion between nevi 
cells and found very little differences with respect to the case in 
which adhesion is absent (see Fig. 8). 

Rupture of the basal membrane. The results reported in 
Fig. 5e and 5f refer to the case in which the basal membrane 
breaks due to the action of mechanical forces. Simulations of the 
rupture of the basal membrane induced by MMPs are reported in 
Fig. 9 and in Movie S8. In particular, the Movie S8 illustrates how 
the effect of MMPs leads to a much rapid invasion of the dermis. 
This is because MMPs typically break the basal membrane in 
different locations, while when the membrane ruptures mechan- 
ically in one location the stress in the rest of the membrane is 
released preventing further rupture. Finally, we notice that the 
average compressive stress in the nevus is reduced when MMPs 
are present (see Fig. 9b). 

Stress dependent growth. Inspired by the experimental 
results on primary melanoma that indicating that compressive 
stress inhibits growth, we perform the simulations with a stress 
dependent division rate, so that compressed cells divide less. The 
morphology of the resulting nevi changes only slightly with respect 
of the previous case (compare Fig. SI with Fig. 5), but we observe 
the development of incipient finger-like structures (see in particular 
Fig.Slc and S If). We also measure the average compressive stress 




100 200 

n 



Figure 6. Mechanical stresses in randomly growing nevi. We report the evolution of the average compressive stress experienced by the 
melanocytes composing the nevus as a function of the size of the nevus, quantified by the total number of cells n, for the same conditions as in Fig. 5. 
The error bars represent the standard error of the mean. The stress is averaged over all the cells in a nevus and over at least ten statistically 
independent realizations of the growth process. The different panels represent different initial locations: a) in the middle of the epidermis, b) in the 
minima of a strong basal membrane, c) in the maxima of a strong basal membrane, d) in the dermis, e) in the minima of a weak basal membrane, f) in 
the maxima of a weak basal membrane. 
doi:1 0.1 371 /journal. pone.0094229.g006 
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Figure 7. Compressive stress distribution in nevi. We report the distribution p(p) of the compressive stresses a (in kPa) experienced by 
melanocytes in a nevus. The results are sampled over different realizations of the process. Since the distribution also changes with time, we only 
consider the time at which the average compressive stress <o-> is highest. The different panels represent different initial conditions corresponding to 
Fig. 5: a) in the middle of the epidermis, b) in the minima of the basal membrane (for strong and weak membranes) c) in the dermis, d) in the maxima 
of the basal membrane (for strong and weak membranes). 
doi:1 0.1 371 /journal. pone.0094229.g007 




250 



Figure 8. Role of adhesion forces and direct interactions with 
the ECM. Morphology of a nevus grown in the dermis a) in presence of 
adhesive forces between nevi cells and b) considering direct 
interactions with the ECM modelled as discussed in the methods 
section, c) The evolution of the average compressive stress experienced 
by the melanocytes composing the nevus for cases a) and b) is 
compared with a simulation in which no adhesive forces and 
interactions with the ECM are present. 
doi:1 0.1 371 /journal. pone.0094229.g008 
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Figure 9. Model of MMP induced breaking of the basal 
membrane, a) We report the typical morphology of a nevus starting 
from the epidermis with a strong basal membrane that can however be 
broken by MMPs. b) The evolution of the average compressive stress in 
presence of MMPs is compared with the case in which MMPs are not 
present. 

doi:10.1371/journal.pone.0094229.g009 
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and find that it is typically smaller than in the previous case (see 
Fig. S2). This is due to the fact that in this case growth occurs in 
regions that are less compressed so that there is more space to 
accommodate new cells. 

Comparison with histological images. In Fig. 10, we 
report two examples of common melanocytic nevi as revealed by 
histological sections. In Fig. 10a, we suggest that an interaepithelial 
juction nevus that is all confined in the epidermis (compare with 
Fig. 5b). Views at different magnifications show that nevi cells are 
closely packed and press against the basal membrane as indicated 
by rounded pocket shown in the left panel. In Fig. 1 Ob, we report 
an example of dermal nevus (compare with Fig. 5d). In this case, 
cells are less packed and are scattered all through the dermis 
without pressing on the basal membrane. These images are in 
good qualitative with the results of our model that show a tendency 
for nevi cells to press from the epidermis into the dermis rather 
than the reverse. This effect is likely due to the different 
mechanical properties of the two skin layers, as assumed in the 
model. 

Discussion 

Melanoma is a very aggressive tumor, being chemio- and radio- 
resistant once it becomes metastatic [44]. Its incidence is 
increasing both in Europe and in US, making it a current 
challenging problem for research. It is commonly believed that 
nevi may progress into melanoma, or alternatively regress by 
differentiating [45]. Since diagnosis is mostly based on dermo- 
scopic and histological analysis, it is extremely important to 
correlate the morphological features of nevi to the degree of 
malignancy or pre-malignancy of the lesions. Since nevi proliferate 
inside a complex microenvironment, we investigated the role of 
mechanics and geometry on the morphology and internal stresses 
of melanocytic nevi. Using human melanoma cell lines, we 
confirm also for melanoma a negative effect of the mechanical 
stress on cellular proliferation as already reported for other tumors 
[8,9,11,12]. We show that osmotic pressure is more effective in 
primary melanoma cells than in metastatic ones, suggesting that 
metastatic cells come from a subpopulation more aggressive and 



insensitive to the mechanical properties of the environment. We 
corroborate this result by experiments of cell growth in collagen 
coated plates where again primary cells are more affected than 
metastatic ones. This suggests that cells composing melanocytic 
nevi may be sensitive to stress as well, an issue that we investigated 
computationally. 

We have devised a computational model for the growth of 
melanocytic nevi in a layered tissue representing the skin. 
Numerical simulations of the model show that the morphology 
of the resulting nevus depends considerably on the environment in 
which the growth takes place, in the epidermis close or far from 
the basal membrane or in the dermis. This is interesting since 
environmental properties, in this case the mechanical behavior of 
the tissue, is rarely considered in the progression of nevi. 
Simulations indicate that nevi are subject to fluctuating compres- 
sive stresses due to the tissue elasticity. If we introduce a stress 
dependence proliferation rate, we observe a decrease in the overall 
stress. Our computational model is appropriate to describe a 
benign nevus, since we do not consider active motion of the cells 
but only their quasi-equilibrium displacement in response to elastic 
forces. To model the growth of melanoma cells, one should 
implement their mobility and their ability to degrade the 
surrounding extracellular matrix, by expressing MMPs. We have 
incorporated this last aspect in the simulations showing that the 
invasion of the dermis occurs more rapidly. It is remarkable that 
many intriguing features observed in melanoma, such as the 
formation of rough tumor boundaries and the invasion of the 
dermis from the epidermis, are observed from simple rules 
combining mechanics and geometry, although we simplified many 
of the irregularities present in the dermis. As a matter of fact we 
expect that, while the overall shape of the nevus can be influenced 
by the degree of accuracy in the dermis description, the qualitative 
picture remains unchanged. A limitation of our model lies in its 
two-dimensional nature that, while it greatly simplifies the 
numerical calculations, could affect the values of the quantities 
we measure. We expect, however, that the general phenomenol- 
ogy we observe should be unchanged in a more realistic three- 
dimensional situation. 




Figure 10. Histological images of melanocytic nevi. We illustrate the morphology observed by optical microscopy at different magnifications 
on histological sections of two types of nevi. a) Intraepithelial junctional nevi are confined in the epidermis and press against the basal membrane 
forming characteristic rounded pockets of densely packed cells (compare with Fig. 5b) b) Dermal nevi are confined in the dermis (compare with Fig. 
5d). Cells are loosely packed and do not touch the basal membrane. The basal membrane is located at the separation between the upper (epidermis) 
and lower (dermis) skin layers and indicated by an arrow in panel a). Images courtesy of Dr. Claudio Clemente. 
doi:1 0.1 371 /journal. pone.0094229.g010 
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Materials and Methods 

Dextran solution 

A master solution of dextran at 10% (w/v) was formed (Dextran 
from Leuconostoc spp, Fluka) and the diluted to the desired 
concentration with complete medium. Transformation from 
dextran concentration to osmotic pressure was performed 
according to the calibration curve measured in Ref. [46]. 

Cell lines 

Human IGR39 and IGR37 cells were obtained from Deutsche 
Sammlung von Mikroorganismen und Zellkulturen GmbH and 
cultured as previously described [47]. IGR39 was derived from a 
primary amelanotic cutaneous tumor and IGR37 from an inguinal 
lymph node metastasis in the same patient. 

Cell growth with collagen 

The cells were plated on pre-coated collagen type I dishes 
(Sigma) according to the manufacturers instructions. 

Colony growth 

Cells are plated on 6 multiwells, fixed with 3.7% paraformal- 
deide (PFA) for 5 minutes and then stained for 30 min with 0.05% 
crystal violet solution. After two washing with tap water, the plates 
are drained by inversion for a couple of minutes. In order to 
control the merging of different colonies, the experiments are 
performed with different initial cell concentrations as described 
recently by our group [42] . Data analysis of the resulting colonies 
has been performed according to Ref. [42]. 

Sulforhodamine B colorimetric assay for cytotoxicity 
screening 

The sulforhodamine B (SRB) assay is used for cell density 
determination based on the cellular protein content according to 
Ref. [48] 

Apoptosis detection 

The Annexin V-FITC Apoptosis Detection kit by Sigma was 
used to detect apoptotic and necrotic cells by detecting annexin V- 
FITG and propidium idodide by flow cytometry, respectively. 
Briefly the cells were incubated with Annexin V-FITC and 
propidium iodide at room temperature for 1 Ominutes and protect 
by light. Then the cells were immediately analysed by flow 
cytometry (FACSAria flow cytometer (Becton, Dickinson and 
Company, BD, Mountain View, CA). Data were analyzed using 
Flowjo software (Tree Star, Inc., San Carlos, CA). 

Histological analysis of bioptic tissues 

Tissue specimens are immediately fixed in neutral buffered 
formalin, embedded in paraffin, sectioned and subjected to 
histopathological characterization after hematoxillin-eosin stain- 
ing. 

Statistical analysis 

Statistical significance is evaluated according to the Kolmo- 
gorov-Smirnov non-parametric test with p< 10 -2 . 

Model simulations 

The model is simulated through a custom made python code. 
To find mechanical equilibrium for the system, we use the fire 
relaxation method [49]. To speed up the code the relaxation 
routine has been written in C. Graphics rendering is done using 



Povray. Simulation and visualization codes are available at 
https:/ / github.com/ alexalemi/ cancersim. 

Supporting Information 

Figure SI Morphology of nevi for different locations of 
the initiating cell for pressure dependent growth. 

Illustration of the results of numerical simulations for nevi grown 
from melanocytes located in different positions in the skin and for 
different mechanical properties of the basal membrane. The 
conditions are the same as in Fig. 5, but here the growth depends 
on the compressive stress acting on each cell. Growing 
melanocytes are shown with a varying color that reflects their 
compressive stresses according to the color bar. The maximum 
<?max and minimum cr m i n for each configuration are equal to (in 
kPa): a) 7.46 10" 6 — 0.13, b) 8.2410" 6 -0.52, c) 0.0016-0.63 d) 
0.00015-0.136 e) 0.00039-0.043 f) 0.0022-0.428. 
(TIFF) 

Figure S2 Mechanical stresses in pressure dependent 
growth of nevi. We report the evolution of the average 
compressive stress experienced by the melanocytes composing 
the nevus as a function of the size of the nevus, quantified by the 
number of cells n, for the same conditions as in Fig. S 1 . The error 
bars represent the standard error of the mean. The different panels 
represent different initial locations: a) in the middle of the 
epidermis, b) in the minima of a strong basal membrane, c) in the 
maxima of a strong basal membrane, d) in the dermis, e) in the 
minima of a weak basal membrane, f) in the maxima of a weak 
basal membrane. 
(TIFF) 

Movie SI Growth of a nevus starting from a melanocyte 
placed deep inside the epidermis. 

(AVI) 

Movie S2 Growth of a nevus starting from a melanocyte 
placed in the epidermis at the minimum of a weak basal 
membrane. 

(AVI) 

Movie S3 Growth of a nevus starting from a melanocyte 
placed in the epidermis at the minimum of a strong 
basal membrane. 

(AVI) 

Movie S4 Growth of a nevus starting from a melanocyte 
placed in the epidermis at the maximum of a weak basal 
membrane. 

(AVI) 

Movie S5 Growth of a nevus starting from a melanocyte 
placed in the epidermis at the maximum of a strong 
basal membrane. 

(AVI) 

Movie S6 Growth of a nevus starting from a melanocyte 
placed deep inside the dermis. 

(AVI) 

Movie S7 Growth of a nevus starting from a melanocyte 
placed deep inside the dermis considering a direct 
interaction with the ECM. 

(AVI) 

Movie S8 Growth of a nevus starting from a melanocyte 
placed in the epidermis at the maximum of a strong 
basal membrane that can be broken by MMPs. 

(AVI) 
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